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The  relationship  between  two  different  descriptions  of  eigenmodes  in 
a  torus  is  investigated.  In  one  the  eigenmodes  are  similar  to  Fourier 
modes  in  a  cylinder  and  are  highly  localised  near  a  particular  rational 
surface.  In  the  other  they  are  the  so-oalied^ ballooning  modes  which 
extend  over  many  rational  surfaces.  Uhm  a  model  which  represents  both 
drift  waves  and  resistive  interchanges  we  investigate  the  transition  from 
one  of  these  structures  to  the  other.  '  In  this  simplified  model  the 
transition  depends  on  a  single  parameter  which  embodies  the  competition 
between  toroidal  coupling  of  Fourier  modes  (which  enhances  ballooning)  and 
variation  in  frequency  of  Fourier  modes  from  one  rational  surface  to 
another  (which  diminishes  ballooning).  As  the  coupling  is  increased  each 
Fourier  mode  acquires  a  sideband  on  an  adjacent  rational  surface  and  these 
sidebands  then  expand  across  the  radius  to  form  the  extended  mode 
described  by  the  conventional  ballooning  mode  approximation.  This 
analysis  shows  tha*  the  ballooning  approximation  is  appropriate  for  drift 
waves  in  a  tokaraafc  but  not  for  resistive  interchanges  in  a  pinch.  In  the 
latter  the  conventional  ballooning  effect  is  negligible  but  they  may 
nevertheless  show  a  ballooning  feature.  This  is  localised  near  the  same 
rational  surface  as  the  primary  Fourier  mode  and  so  does  not  lead  to  a 
radially  extended  structure.  • — - * — — • — - tS-""' 
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1.  INTRODUCTION 


Two  essentially  distinct  descriptions  have  been  used  for  short  wave 
length  perturbations,  such  as  drift  waves  and  resistive  instabilities,  in 
an  axisymmetric  toroidal  plasma.  The  first  regards  toroidal  curvature  as 
a  small  effect  and  expresses  the  plasma  perturbation  in  terms  of  coupled 
cylindrical  Fourier  modes  with  toroidal  and  poloidal  mode  number  n  and 
m.  In  the  absence  of  toroidal  effects  each  Fourier  mode  is  centered  on 
its  appropriate  rational  surface  where  m  =  nq(r)  and  its  properties  are 
determined  by  the  magnetic  shear,  density  gradient  etc.  at  the  rational 
surface;  derivatives  of  the  shear  and  density  gradient  are  neglected.  In 
this  approximation  modes  with  different  m  are  almost  degenerate.  In  the 
case  of  drift  waves  an  "outgoing  wave”  boundary  condition  applied  to  each 
Fourier  mode  leads  to  the  so-called  'shear  damping'  effect.1 

On  the  other  hand,  in  the  second  description  it  is  precisely  the 

degeneracy  of  the  Fourier  modes  which  is  exploited.  This  leads  to  the 

2 

ballooning  representation  in  which  disturbances  centered  on  neighboring 
rational  surfaces  are  strongly  coupled  to  produce  a  radially  extended 
structure  -  similar  to  the  quasi-mode  of  Roberts  and  Taylor3  -  which 

"balloons"  on  one  side  of  the  torus.  The  coupling  reduces  the  shear 

4 

damping  of  drift  waves  -  and  it  may  disappear  altogether  if  the  coupling 
is  sufficiently  strong. 

In  this  paper  we  discuss  the  relationship  between  these  two  opposing 
views  of  mode  structure  in  toroidal  systems.  In  Section  2  we  introduce  a 
simple  model  for  drift  waves  in  which  Fourier  modes  are  coupled  but  not 
degenerate.  In  Section  3  we  show  that  a  similar  model  can  describe 


resistive  interchange  instabilities.  In  Sections  4  and  5  we  use  these 
models  to  explore  the  transition  from  a  Fourier  mode  description  to  a 
ballooning  mode  description  as  the  toroidal  curvature  is  increased. 

This  investigation  shows  that  as  the  transition  occurs,  each 
Fourier  mode  first  acquires  a  ballooning  component  centered  on 
neighboring  rational  surfaces  and  this  structure  then  expands  across  the 
radius  to  form  the  conventional  ballooning  mode.  However,  for  resistive 
interchanges  the  overlap  between  Fourier  modes  on  adjacent  rational 
surfaces  is  exponentially  small  and  this  extended  ballooning  structure 
does  not  arise.  Instead  an  alternative  ballooning  effect,  discussed  in 
Section  6,  can  occur.  In  this  the  ballooning  contribution  is  centered  on 
the  same  rational  surface  as  the  primary  Fourier  mode.  Consequently 
although  the  mode  "balloons"  it  does  not  have  the  quasi-mode  character. 

In  Section  7  we  discuss  the  relationship  of  these  results  to  the 
conventional  ballooning  mode  approximation. 


2 .  DRIFT  WAVES 


A  simple  model  for  drift  waves  in  a  large  aspect  ratio  tokamak  has 
been  described  in  Refs.  5,  6  and  7.  The  perturbed  potential  is  written  in 
the  form: 


$  =  <j>(  0.  p)  exp  ( i < n C  -  nigS)  -  iwt  ]  (1) 

where  p  is  the  radial  distance  from  the  rational  surface  m^  =  nq(r) 
and  8  is  the  poloidal  angle.  Then  $( 9, p)  satisfies 
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<Lq  +  L1  +  L2 )  <t>  =  0 


(2) 


where 


L  =  a?  *1  -  b 

1  dp2 


{ oi  —  ui.  -  ioi6)  ai.  e  o  . 


-  (—  •  _)  (_  +  iksp)  ,  (3) 
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1  *  r  „  .  isin0  6  > 

2  e  _  cos  9  + _ 

n  v  a  ' 

to  k  op 
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L2  (arc  +  a), ) 


(5) 


In  Eqs.  13)  —  (5),  k  =  nq/r  ,  s  =  (r/q)(dq/dr)  ,  a2  =  c^itTVe^2  , 
b  «  k2a2,  en  =  qec  =  r^/Rt  “here  is  the  density  scale  length 

and  t  *  Te/Ti  ■  The  diamagnetic  frequency  a).  =  (kcTe/eBrn)  is  expanded 
about  a  local  maximum  with  L~2  =  (l/a>t )  (d  2a)t/dr 2 )  and  the  parameter  i6 
represents  the  destabilising  effect  of  electron  Landau  resonance  and 

trapped  electrons.  This  model  can  be  obtained  from  the  general  equation 

0 

for  drift  waves  in  a  tokamak,  given  by  Tang,  by  taking  the  long 
wavelength  limit  ka^  <<  1  and  neglecting  temperature  gradients. 

The  operator  LQ  describes  drift  waves  in  a  cylinder  (e.g.  as 
discussed  by  Pearlstein  &  Berk^)  and  depends  only  on  the  shear  and  density 
gradient  at  the  rational  surface  itself;  L^  describes  the  toroidal 
effects  and  L^  introduces  the  radial  variation  of  . 

To  simplify  the  problem  further  we  consider  the  regime  ecs  , 
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b  ,  6  <  1  ,  in  which  oo  =  .  Then  Eq.  (2)  becomes 


*2  »■  ft  ft 

[ —  -  0*[ —  +  ix)2  -  e(cos6  +  is  sin6  — )  -  <x2  -  =  0  (6) 

3x2  39  5x 

where  x  =  kps  .  The  three  parameters  in  the  model  are  then 

e  2e 

c  n  1 

o  =  —  ,  e  =  -  ,  k  =  _  (7) 

bs  bs2  k^^bs^d  +  t) 

and  the  eigenvalue  \  is  related  to  the  mode  frequency  u  by 

\  =  - ] -  tiL  -  1  +  b(1  +  t)  -  i6)  .  (8) 

bs  2(  1  +  t)  * 

In  the  approximation  w  =  a the  parameter  a 2  is  real  and  there  is 
no  explicit  dissipation  in  the  model.  However,  in  the  usual  way,  we 
suppose  that  to  has  a  small  positive  imaginary  part  in  order  to  justify 
neglect  of  initial  value  terms.  Correspondingly  a2  acquires  a  small 
negative  imaginary  part. 


3.  RESISTIVE  INSTABILITIES 


Another  instability  which  can  be  represented  by  a  model  similar  to 

3  9 

Eq.  (6)  is  the  localised  electrostatic  resistive  interchange  '  in  a  large 
aspect  ratio  toroidal  pinch.  In  this  case  the  model  can  be  derived  from 
the  high  mode  number  resistive  ballooning  equation  discussed  by  Bateman 
and  Nelson.10  In  the  electrostatic  limit  y  <  gn2|Vg|2  [where  Y  is  the 
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growth  rate,  T)  the  resistivity  and  the  mode  has  an  eikonal 
representation  $  ~  p  exp(inS)  ]  this  equation  takes  the  form 

|  VS !  2  B  x  vs  *K 

_I_  B*V  _  B •  V4>  -  p„v2  -  0  +  2p'  2 - 2  p  =  0  (9) 

TP2  ~  B2  ~  °  B2  B2 

where  is  the  density  and  K  is  the  field  line  curvature. 

In  a  large  aspect  ratio  pinch,  with  magnetic  field 
B  =  BQ(r)[l  -  A(r)  cosQ]  ,  Eq.  (9)  can  be  further  simplified. 

Reverting  to  the  form  ( 1 )  for  the  perturbations  and  taking  the  long 
wavelength  limit  it  becomes 

[— —  +  <j2(—  +  ix  )  +  a  (p  cos  6  +  isp  sin0-5_)  -Kx2+a]i)>  »  0 

dx2  R  36  1  2  3x  R 

(10) 
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and  L~  2  =  -  ( 1/a) (d 2at/dr 2)  represents  the  radial  variation  of  the 
cylindrical  curvature  term  a  . 


4.  DEGENERACY  AND  ITS  REMOVAL 


We  have  noted  in  Sections  2  and  3  that  drift  waves  and  resistive 
instabilities  in  a  large  aspect  ratio  torus  can  both  be  described  by  the 
model  equation 


[— —  -  a2  (^_  +  ix)2  -  e(cos9  +  is  sin6  — )  -  <x2  -  \J<p 
dx 2  d0  5x 


0 


(12) 


(where  we  set  ap1  =  ap2  =  e  in  the  resistive  case).  If  Fourier  modes 
~  exp(-im6)  are  introduced  into  Eq.  (12)  it  can  be  seen  that  a 
determines  their  radial  width,  k  describes  the  difference  between  modes 
on  neighboring  rational  surfaces  and  e  describes  the  coupling  between 
modes.  Thus,  despite  the  many  simplifications,  the  model  retains  the 
features  necessary  for  both  Fourier  mode  and  ballooning  mode  descriptions. 
The  distinction  between  drift  and  resistive  instabilities  lies  in  the  sign 
of  a2  .  For  drift  waves  a2  >  0  and  for  resistive  instabilities 
a2  <  0  .  In  this  section  we  concentrate  on  drift  waves. 


We  treat  both  e  and  <  as  first  order  small  quantities  and  seek  a 
solution  of  Eq.  (12)  in  the  form 
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where 


u^(x)  =  u(x  -  m)  =  exp[-io(x  -  m)2/2]  (14) 


satisfies  the  zero-order  equation 


with  \Q  =  -  io  . 

Each  u  (x)  represents  a  separate  wave  propagating  from  its  rational 
m 

surface  x  =  m  .  When  a  has  the  small  imaginary  part  referred  to  in 
Section  2  the  amplitude  of  this  wave  decreases  with  distance  from  the 
rational  surface.  The  eigenvalue  \ corresponds  to  a  drift  wave 
frequency 


b(1  ♦  x)  +  i[6  -  C-~)] 


which  exhibits  the  competition  between  the  destabilising  term  6  and  the 

•shear  damping'  ~  (r  s/Rq)  associated  with  the  'outgoing  wave' 
n 

feature^)  of  the  u  (x)  . 

m 

So  far,  the  modes  u  (x)  are  degenerate  and  the  coefficients  c 
m  m 


are  completely  arbitrary.  The  appropriate  combination  of  c  for  the 

m 

toroidal  problem  is  determined  when  the  additional  terms  of  Eq.  (12)  are 
introduced  and  the  degeneracy  is  removed.  Then  vre  have 


Vi  "  l€L1  +  KL2>*0  "  Vi  "  Vo 


(17) 


« 

where  LQ  is  the  differential  operator  which  appears  in  Eq.  (15)  and  el^ 

A 

and  kL^  are  the  remaining  operators  of  Eq.  (12).  The  quantity  can 

be  annihilated  by  the  operator 


/  dx  4  d8  u  (x)  exp(imG) 
'  '  m 


(18) 


This  leads  to  a  recurrence  relation  for  the  c  : 

m 


(\,<u2>  +  K  <x2u2>  )  c  +  (<u  u  J ,>  +  s  <u  Du  >  )  c 

1  m  m  >  m  2  *■  m  m+1  m  ra+1  '  m+1 

+  4-(<uu  -  s  <u  Du  .>]  c  =  0 

2  v  m  m-1  m  m-1  '  m-1 


(19) 


where 


<...>  =  /  ...  dx  and  D  =  (20) 

Ox 


Using  the  explicit  expression  (14)  for  the  u  (x)  ,  Eq.  (19)  can  be 

m 

reduced  to  a  form  well-known  in  the  theory  of  Mathieu  functions: 
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(21  ) 


Note  that  the  symbol  q  is  used  here  to  conform  with  standard  notation 

for  Mathieu  functions:  it  should  not  be  confused  with  the  toroidal 

'safety  factor'  in  Section  1.  The  boundary  conditions  on  Eq.  (21)  are 

that  c  ■¥  0  as  I  ml  +  ®  and  the  required  solutions  are  related  to  the 
m  I  1 

' period- u'  Mathieu  functions  ce  and  se 


S.  STRUCTURE  OF  TOROIDAL  EIGENMOPES 

Tne  preceding  section  shows  that  the  structure  of  the  toroidal 
eigenmodes  is  determined,  through  the  recurrence  relation  (21),  by  the 
parameter  q  .  This  parameter  incorporates  the  influence  of  toroidal 
curvature,  the  width  of  a  basic  Fourier  mode  and  the  variation  of 
diamagnetic  frequency  with  the  location  of  its  resonant  surface. 

When  jqj  is  small  the  fundamental  eigenvalue  of  Eq.  (21)  is11 

a  =  -2q2  +■  1  q1*  +  ...  (23) 

and  the  associated  coefficients  c  decrease  rapidly  with  m  , 

m 
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'0 


1 


-2q  , 


(24) 


Consequently,  when  |q|  is  small  only  is  significant.  Recalling 

that  a  factor  ~  expl-im^S)  was  extracted  from  the  perturbation  (Eq.  (1)} 
this  shows  that  the  toroidal  eigenmode  comprises  a  single,  well-localised 
Fourier  mode  on  the  rational  surface  nq{r)  =  m^  with  weak  and  rapidly 
decreasing  sidebands  on  adjacent  rational  surfaces.  These  weak  sidebands 
provide  the  only  ballooning  effect.  There  is  a  small  reduction  in  shear 
damping  (compared  to  Eq.  ( 16)  )  but  this  is  due  to  the  variation  in 
across  the  plasma  profile  rather  than  to  toroidal  coupling. 


On  the  other  hand,  when  |q|  is  large  the  toroidal  eigenmodes  take 
on  an  entirely  different  structure.  In  one  of  them  c^  varies  slowly 
with  m  and  Eq.  (2i)  may  be  replaced  by  a  differential  equation,  treating 

(4) 

m  as  a  continuous  variable  : 

d2c 

q  - -  +  (m2  -  a  +  2q)  c  =  0  (25) 

.  2  m 
dm* 

Then  the  fundamental  eigenvalue  is 

a  =  2q  -  iq1/2  (26) 
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from  m  =  0  out  to  m  ~  ±jq|  and  fall  off  exponentially  beyond  this 


point. 

There  is  also  another  eigenmode  when  q  is  large.  In  this  c  ~ 

m 

exp {iron)  -  which  corresponds  to  replacing  the  poloidal  angle  0  by 

0+7i  in  a  mode  with  slowly  varying  c  .  This  signifies  that  this 

m 

second  mode  is  centered  on  the  inside  edge  of  the  torus  (  0  =  it)  rather 

than  the  outside  (0  »  0).  Writing  c  »  elmilc'  ,  where  c'  varies 

mm  m 

slowly,  one  finds  for  this  eigenmode 

a  =  -2q  +  q1/2  (28) 


and 


=  exp(imtt  -  m2/2q1/2)  (29) 

with  -n  <  Arg  q  <  n  . 

At  large  q  ,  therefore,  the  toroidal  eigenmodes  extend  over  many 
rational  surfaces,  involve  many  poloidal  m  numbers  and  have  a  strong 
ballooning  character .  This  is  the  situation  described  by  the  conventional 
ballooning  mode  approximation.  If  Re(q)  >  0  the  more  strongly 
ballooning  mode  is  centered  at  9  =  Tt  ,  (on  the  inside  edge  of  the  torus) 
and  the  shear  damping  of  this  mode  is  reduced  by  toroidal  coupling  if 
s  >  1/2  -  in  agreement  with  the  detailed  computations  of  toroidal  drift 
modes  by  Hastie,  Hesketh  and  Taylor.7 


-1  I- 


In  a  typical  tokamak  {  B  ~  3T  ,  ~  T\  ~  2  keV  #  a  ~  20  cm  , 

R  ^  100  cm  )  q  is  ~  10  ^  -  10  ^  for  drift  waves  with  ka .  ~  1  so  that 

l 

such  drift  waves  lie  well  within  the  ballooning  regime.  They  are 
therefore  well-described  by  the  conventional  ballooning  mode  approximation 
(see  Section  7). 


6.  ANOTHER  BALLOONING  EFFECT 


For  resistive  interchanges  we  replace  o  by  -iaR  and  t^ie  coupling 
coefficient  q  becomes 


€R  {  °RS  \  ,  aR  x 

(1  f  -j-)  exp(-  _) 


(30) 


where  we  have  again  set  p1  =  and  eR  =  ap^  . 

In  a  typical  pinch  experiment  aR  is  large  and  qR  is  exponentially 
small.  Consequently  resistive  interchanges  do  not  show  the  ballooning 
effect  discussed  so  far  and  are  not  described  by  the  conventional 
ballooning  mode  approximation.  However,  they  may  show  a  different 
ballooning  effect  which,  although  small  when  cjR  is  large,  is  not 
exponentially  small. 


To  investigate  this  alternative  ballooning  effect  we  return  to  the 
model  equation  (10) 
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[-5 —  +  a2  (-5-  +  ix )  +  a(l  +  p1  cos6  +  isp2  sing  -xx2] 


where  we  drop  the  subscripts  on  a  and  x  but  explicitly  indicate  the 
normal  and  geodesic  curvatures  since  the  geodesic  curvature  p^  plays  a 
dominant  role  when  a  >>  1  • 


We  now  observe  that  when  a  >>  1  the  basic  Fourier  mode  is  very 
narrow  and  therefore  introduce  a  new  length  scale  y  =  o1/2x  ;  then 


ae2  o1'2 


i  ,  a 
~  (2y  <55  +  — 


+  "5  *  y2^  +  +  -5  +  P-, cos Q ) 3 1 


We  now  seek  a  solution 


ly.9)  +  -  $..(y,6)  +  ... 

0  r.l/2  1 


and,  to  ensure  maximal  ordering,  we  take  a  ~  a  ,  x  ~  o2  .  Then  in 

lowest  order  4  =  4  (y)  and  in  first  order 

0  0 


iosp2  d$Q 

__  Sin0^—  +  $l(y)  . 


In  second  order,  after  annihilating  by  integration  over  0  ,  we  have 


3y  ^ 


o 


(35) 


Thus  at  large  o  ,  differs  slightly  from  the  basic  Fourier  mode 

and  there  is  a  small  shift  in  the  eigenvalue.  However,  the  interesting 
feature  is  the  appearance  of  the  contribution  ,  proportional  to  sin6  . 
This  is  a  "ballooning"  component  located  near  the  same  rational  surface  as 
the  primary  Fourier  mode  4^  and  proportional  to  cr322  .  it  is  thus 
quite  distinct  from  the  ballooning  components  considered  hitherto  which 
are  located  on  rational  surfaces  adjacent  to  that  of  the  primary  mode  and 
are  exponentially  small  in  a  .  It  is,  however,  the  effect  referred  to  as 
"ballooning"  in  discussions  of  interchange  modes  by  some  Soviet  authors 
(see,  for  example,  Ref.  12). 


This  'centered'  ballooning  effect  can  be  regarded  as  a  result  of 

coupling  between  the  harmonics  of  the  basic  Fourier  modes  u  which  are 
-  m 

usually  neglected.  These  harmonics  are  given  by  Hermite  functions 

uk(x)  *  (_j5 - )  exp [-  5.  (x  -  m)2]  H  f/o(x  -  m)  ]  (36) 

”  2*kl/S  *  K 

with  k  >  0  .  If  we  return  to  Section  4  and  include  these  harmonics  in 
the  expansion  of  4^  ,  so  that 

k  k 

”  I  l  c  u  <x)  exp(-im9)  (37) 

0  ,  m  m 

k  m 

then  in  place  of  the  recurrence  relation  (21)  we  obtain  (retaining  only 
the  dominant  geodesic  curvature  in  Eq.  (31)  and  ignoring  x  ), 
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,  .  k 

(\  '  0  °m 


Compared  to  the  earlier  calculation  of  Section  4,  there  is  a  change  in  the 


coefficient  of 
amount 


equivalent  to  a  shift  in  the  eigenvalue 
m 


a  by  an 


‘S2 


1*0 


( V^o1 


{<u°  D  U 
1  m  m+l 


+  <uw  D  u 


m- 1 


(42) 


Eq.  (40)  shows  that  there  is  also  an  additional  component  of  the 
perturbation,  given  by 


VP2 

2 


-i(m+1 ) 0 


1 

1*0 


1  „  i 

, - r  <u 

Ta^-aT  m+1 


D  uu> 
m 


m+ 1 


(43) 


The  individual  matrix  elements  in  (42)  and  (43)  involve  the  overlap 
between  Hermite  functions  centered  on  different  rational  surfaces  and  are 
consequently  exponentially  small  when  a  >>  1  •  However,  the  expressions 
(42)  and  (43)  are  not  themselves  exponentially  small.  It  can  be  shown, 
using  the  Hermite  generating  functions,  that  for  large  o  the  expression 
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(42)  is 


“o®  2p2//  °2  and  agrees  with  the  shift  in  eigenvalue  given  by  Eq. 

(35).  Similarly  the  expression  (43)  is  ~  OgSp^a  and  agrees  with  Eq. 

(34)  at  large  a  .  Note  again,  that  although  this  contribution  to  the 

perturbation  is  proportional  to  exp (-i(m+1 ) 6)  it  is  localised  near  the 
th 

m  rational  surface,  (i.e.  near  the  same  surface  as  the  primary 

component  c  0  ) . 

m 

7.  THE  BALLOONING  REPRESENTATION 

2  7 

In  the  usual  theory  of  high  n  ballooning  modes  '  one  maps  the 
poloidal  angle  9  on  to  an  extended  coordinate  ti  with  -»  <  ri  <  »  an(j 
writes  the  perturbation  in  the  form 

<j>(x,  n>  =  Mx)  e  ix(n  +  k)  .  (44) 

Then  A(x)  is  taken  to  be  a  slowly  varying  ’envelope’  (on  the  scale 
xTii  )  and  to  lowest  order  in  (1/i/n  >»  f(T),x)  satisfies  a  ’local’ 

one-dimensional  eigenvalue  equation  in  the  extended  q  coordinate  (with 
x  fixed).  For  the  present  drift  wave  model  this  one-dimensional  equation 
is 

[a2  — —  +  ( h  +  k)  2  +  e(cost)  +  s(  ri  +  kjsinri)-  kx  2  +  \Jf  ( t),x)  »  0  ,  (45) 

dri2 
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where  X(k,x)  is  the  local  eigenvalue. 


The  parameter  k  ,  representing  an  as  yet  undetermined  radial 
wave-number,  and  an  equation  for  the  radial  envelope  function  A(x)  are 
determined  in  higher  orders  of  the  expansion  in  ^^//n)  .  Thus  k  is 
obtained  from  the  equation 


3\ 

8k 


(x,k) 


0 


(46) 


and  A(x)  satisfies  the  equation 


1  a2X  d2A  ,  x2  a2xu 

- +  (\  -  x  - - Ja 

2  3k2  dx2  °  2  dx2 


(47) 


The  condition  for  self-consistency  of  the  theory,  that  A(x)  vary 
slowly,  is  satisfied  therefore  if 


8k2  dx2 


(48) 


To  assess  this  criterion  in  the  present  model  we  can  determine 
X(k,x)  from  Eq.  (45)  by  perturbation  in  e  and  <  .  [This  expansion  is 
quite  distinct  from  the  ballooning  expansion  in  1/i^n  .  ]  In  leading  order 
the  outgoing  wave  condition7  implies 


fQ  =  exp [ i ( r)  +  k)  2/2o] 
with  eigenvalue  XQ  =  -io  .  In  next  order 


(49) 
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[°2  —  +  <n  +  k) 2  +  ^  ]f, 
dri2 

■ 

+  [e(cosr)  +  s(r)  +  k)sinT))-  icx2  +  ]f  =  u 

CO 

Annihilating  by  the  operation  j  drjfg  (recalling  that  a  nas  a 

—  00 

small  positive  imaginary  part)  we  find 

\  =  kx2  -  2<q  cosk  (51, 

where  q  is  the  parameter  introduced  in  Section  4. 

The  requirement  of  Eq.  (46)  can  therefore  be  satisfied  by  k  =  0  or 
k  =  re  ,  corresponding  to  ballooning  modes  centered  on  the  inside  or  the 
outside  of  the  torus,  consistent  with  the  discussion  of  Section  5.  It 
then  follows  from  Eq.  (51)  that 

— -  =  t  2kq  and  -2.  _  =  2  <  (52) 

t)k  2  dx  2 

and  the  self-consistency  condition  (48)  reduces  to  |q|  >>  1  •  This  is 
again  precisely  what  we  expect  in  the  light  of  the  discussion  given  in 
Section  5  and  confirms  that  the  ballooning  approximation  is  valid  when 

H  »  1  • 

When  |q|  is  small  the  ballooning  approximation  is  invalid;  A(x) 
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does  not  vary  slowly  and  the  consistency  condition  (48)  is  not  satisfied. 
Nevertheless  the  ballooning  formalism  may  sometimes  be  useful  even  when 
|q|  <  1  .  This  is  because  the  Ansatz 


<J>( x,  t))  =  e  iXT1f  (  T)) 


(53) 


> 


satisfies  the  model  equation  exactly  if  ic/e  can  be  neglected.  For 
drift  waves  this  is  never  possible  when  |q|  <  ^  since  </ e  must  then 
be  large.  However,  for  resistive  interchanges,  when  q  is  given  by  Eq. 
(30)  and  aR  is  large,  we  may  have  both  |q|  and  k/ c  small.  In  such 
an  event  the  lowest  order  approximation  of  ballooning  theory  (analogous  to 
Eq.  (45))  provides  a  good  estimate  for  the  correct  eigenvalue  (and  f(r|) 
incorporates  the  alternative  ballooning  effect  discussed  in  Section  6)  - 
even  though  the  formal  development  of  the  ballooning  approximation  breaks 
down. 


8.  SUMMARY  AND  CONCLUSIONS 


We  have  investigated  the  connection  between  two  seemingly  opposed 

views  of  toroidal  eigenmodes.  In  one  the  eigenmodes  are  basically 

cylinder  Fourier  modes  ~u  (x)exp(-im8)  localised  near  a  particular 

in 

rational  surface.  In  the  other,  the  eigenmodes  are  ballooning  structures 
which  extend  over  many  rational  surfaces.  Using  a  model  which  represents 
both  drift  waves  and  electrostatic  resistive  interchange  instabilities,  we 
have  determined  the  parameter  which  determines  whether  the  eigenmodes  are 
Fourier-like  or  ballooning-like .  This  parameter  q  ,  defined  in  Eq.  (22), 
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embodies  the  competing  physical  effects  of  coupling  between  Fourier  modes 
(which  enhances  the  ballooning  tendency)  and  variation  in  the  frequency  of 
Fourier  modes  from  one  rational  surface  to  another  (which  reduces  the 
ballooning  tendency). 

As  q  increases,  the  conventional  ballooning  effect  first  manifests 

itself  through  coupling  of  the  fundamental  Fourier  modes  u°(x)  on 

m 

adjacent  surfaces*  This  produces  a  toroidal  eigenfunction  consisting  of  a 
primary  Fourier  mode  with  weak  'sidebands'  on  neighboring  rational 
surfaces.  These  'sidebands'  then  increase  in  amplitude  and  extend  across 
the  radius  until  the  'quasi-mode ' -like  structure  of  the  conventional 
toroidal  ballooning  mode  is  produced.  When  |q|  is  large  the  ballooning 
mode  extends  over  ~  qA/l*  rational  surfaces. 

2 

This  analysis  supports  the  criterion  given  earlier  for  the  validity 
of  the  conventional  ballooning  approximation  in  a  toroidal  system  and 
shows  that  this  approximation  is  appropriate  for  drift  waves  in  a  tokamak 
but  not  for  resistive  interchanges  in  a  pinch.  For  resistive  interchanges 
the  conventional  ballooning  effect  is  negligible  but  they  may  show  another 
ballooning  effect.  This  involves  a  sideband  centered  on  the  same  rational 
surface  as  the  primary  Fourier  mode  and  does  not  lead  to  a  radially 
extended  quasi-mode  structure.  It  can  be  considered  as  the  result  of 
indirect  coupling  between  harmonics  of  the  Fourier  modes.  Surprisingly, 
although  the  conventional  ballooning  theory  then  breaks  down,  its  lowest 
order  approximation  may  still  be  useful. 
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